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Abstract 

We study the four-dimensional Yang-Mills theory in the presence of a three-dimensional 
membrane of fermions by lattice Monte Carlo simulations. We analyze the phase 
structure of this theory at finite temperature. Below the phase transition tem- 
perature of the pure Yang-Mills theory, we obtain an unconventional phase with 
spatially-nonuniform vacuum. In this phase, the expectation value of the Polyakov 
loop is finite on the membrane, and it exponentially decays to zero outside the 
membrane. 
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1 Introduction 



As a theoretical interest in the gauge field theory, we can consider the theory 
which contains the gauge field and the matter field in different dimensions. 
Such a theory often shows characteristic properties [1,2,3.4]. In condensed mat- 
ter physics, it has been realized in trapped electron systems, such as graphene 
[5]. Recently, graphene is investigated in lattice gauge theory [rJfTfgfl?] . In the 
gauge/string duality, this kind of setup has been frequently used in terms of 
D-brane ITUl. 



Motivated by these works, we consider the Yang-Mills theory coupled with 
fermions. We depict our theoretical setup schematically in Fig. [TJ The Yang- 
Mills field lives in the (3+l)-dimensional space-time, while the fermion field 
lives only on the (2+l)-dimensional hyperplane at z = 0. The existence of 
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Fig. 1. The schematic figure of the (3+l)-dimensional gauge theory with (2+l)-di- 
mensional fermions. The fermions are localized at z = 0. 

this fermion membrane breaks the Lorentz invariance and the translational 
invariance in the ^-direction. From the viewpoint of the gauge field, the fermion 
field is localized in the low-dimensional subsystem. From the viewpoint of the 
fermion field, the background gauge field has one extra dimension. 

This theory is renormalizable. Counter terms are introduced for the gauge field 
in four dimension and for the fermion field in three dimension at z = 0. There 
is no additional ultraviolet divergence, unlike two-dimensional fermions in the 
four-dimensional Yang-Mills theory [2ll4] . This theory is classically conformal 
when the fermions are massless. There is no dimensional parameter except for 
the bare fermion mass, unlike the three-dimensional Yang-Mills theory. The 
dynamics is uniquely determined through dimensional transmutation. 

In this Letter, we study this system by means of lattice Monte Carlo simula- 
tions. We focus on how the fermion membrane affects the vacuum structure 
at finite temperature. In the finite-temperature formalism, the t-direction is 
compactified in the Euclidean metric with a periodic boundary condition for 
the gauge field and with an antiperiodic boundary condition for the fermion 
field. 



2 Formalism 

This system is formulated by the standard technique of lattice gauge theory. 
For the gauge part of the lattice action, we used the Wilson gauge action 




(1) 



2 



where U^ v {x) is the plaquette variable, i.e., the minimal gauge invariant loop. 
The lattice spacing a is determined by the dimensionless parameter = 
2N c /g 2 . Because this system is equivalent to the pure gauge theory in large z, 
we can set the same physical scale as in the pure gauge theory. For the fermion 
part, we used the three-dimensional staggered fermion, 



where x{ x ) * s the spinorless staggered field and rj^x) is the staggered phase 
The summations are taken over the three-dimensional plane at z — 0. In 
the continuum limit, this fermion action corresponds to the standard three- 
dimensional form X) 4 > {o'iD i ±m)ip and the fermion field is the four-flavor two- 
component spinor. The four flavors contain two positive-mass states of +m and 
two negative-mass states of — m. This fermion field is equivalently rewritten as 
the two-flavor four-component spinor under a 4 x 4 representation of the Dirac 
matrices [TTp2] . This four-component spinor construction preserves parity and 
time-reversal invariances. 

To generate the gauge configurations under the full lattice action Sq + Sf, 
we made use of the Hybrid Monte Carlo algorithm. The gauge configuration 
with the lattice volume N% x N t includes the dynamical fermions on the three- 
dimensional volume Ng x N t . The parameters are set at N c = 3 and — 5.7, 
and the corresponding lattice spacing is a ~ 0.19 fm. The bare fermion mass 
is set at ma = 0.2. To simulate several values of temperature, we changed the 
temporal extent Nt with the fixed lattice spacing. The physical temperature 
is given as T = l/(N t a) with a" 1 ~ 1 GeV. The lattice sizes of the gauge 
configurations are listed in Table [TJ 



3 Results 

First, we measured the fermion condensate 



at the z = plane. We show the numerical data in Table [T] and Fig. [2J 
Note that the fermion condensate is always finite because of the finite fermion 
mass. The fermion condensate decreases as temperature increases, and its 
derivative seems to be large in 125 MeV < T < 167 MeV. This suggests that 
the remnant of spontaneous symmetry breaking contributes to the fermion 
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Table 1 

Numerical data of simulations. The temporal extent Nt, the spatial extent N s , the 
physical temperature T, the fermion condensate £ and the Polyakov loop P are 
listed. The statistical errors are shown in parentheses. 
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0.734(12) 


0.0017(22) 


-0.0011(28) 
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0.692(14) 


0.0065(34) 


0.0007(27) 
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0.596(13) 


0.0252(31) 


-0.0051(30) 
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0.1025(44) 


0.0121(41) 
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0.378(8) 


0.1128(46) 


-0.0030(33) 
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0.379(6) 


0.0999(30) 


-0.0008(22) 
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16 


250 


0.228(1) 


0.2966(38) 


0.1538(43) 



condensate in low temperature and it is restored in high temperature. The 
broken symmetry is flavor symmetry, not chiral symmetry, because there is no 
chiral symmetry in three dimension. In the massless and continuum limit, the 
symmetry breaking pattern is considered to be U(Nf) —> U(Nf/2) x U(Nf/2), 
as three-dimensional QCD [T31[T4"] . 

Next, we calculated the expectation value of the Polyakov loop 



which is a good indicator for confinement. Unlike the fermion condensate, 
we can define the Polyakov loop not only on the fermion membrane (z = 0) 
but also outside it (z ^ 0). We show the Polyakov loop value at z = 
and z = aN s /2 in Table [1] and Fig. [2j In large z limit, this theory should 
be equivalent to the pure Yang-Mills theory. The SU(3) pure gauge theory 
has the first-order phase transition at T c ~ 250 MeV, e.g., /3 C ~ 5.69 with 
N t = 4 [15J The data of z = aN s /2 is consistent with this expectation. On 
the other hand, the Polyakov loop at z = shows a different behavior. It is 
finite even at T ~ 125 MeV. As a result, we obtain an unconventional phase in 
125 MeV < T < 250 MeV. In this phase, the Polyakov loop value is nonzero 
at z — 0, whereas it is almost zero in large z. This phase is interpreted as 
confinement phase with partial deconfinement around the fermion membrane. 

The existence of this phase is understood as follows. Apart from the mem- 
brane, the system has the first-order phase transition of the pure Yang-Mills 
theory at T c ~ 250 MeV. On the membrane, however, the theory includes the 
dynamical fermions. This makes the coupling constant smaller and the transi- 
tion temperature lower locally around the membrane. Therefore, the fermion 
membrane induces partial deconfinement below the transition temperature of 
the pure Yang-Mills theory. 
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Fig. 2. The Polyakov loop P and the fermion condensate £ with N s = 16. 
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Fig. 3. The volume dependence of the Polyakov loop P and the chiral condensate 
£ with N t = 6. 

In the present simulations, we cannot determine the temperature and the 
order of the phase transition. These properties would depend on the fermion 
mass and the flavor number. To validate the existence of this phase, we here 
check the dependence on the spatial volume V = a 3 N 3 at fixed temperature. 
In Fig. [31 we plot the fermion condensate and the Polyakov loop at T ~ 167 
MeV as a function of the inverse volume 1/V. All these are insensitive to the 
volume. In particular, the Polyakov loop at z = seems to remain finite in 
the infinite volume limit 1/V — > 0. Thus, at least, the calculation at T ~ 167 
MeV, i.e., N t = 6, lies in this phase. 



In this phase, the vacuum structure is spatially-nonuniform in the z-direction. 
In Fig. U we show the Polyakov loop value as a function of z. Since the 
distribution is symmetric about z = 0, we only show the region of < z < 
aN s /2. The numerical data is well fitted by Cexp(— z/zq) + Pq. At T ~ 167 
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Fig. 4. The z-dependence of the Polyakov loop P with N s 
best-fit functions of Cexp(— z/zq) + Pq. 



16. The curves are the 



MeV, i.e., N t = 6, the Polyakov loop value is nonzero at z = and decreases 
exponentially in z > 0. The best-fit damping parameter is zq ~ 0.3 fm. Taking 
into account the region of z < 0, the thickness of the deconfinement layer is 
roughly 1 fm. 

As another observable, we consider the color averaged potential between static 
color sources. The color averaged potential V(R) is extracted from the corre- 
lation function of the Polyakov loop and the anti-Polyakov loop, 



v(R) 

e T 



t-2 



(5) 



with R — \xi~ x 2 \. In the case of N c = 3, the color averaged potential includes 
the color-singlet and color-octet components. In a deconfinement medium, the 
potential is Debye screened as the color sources are separated. We measured 
the color averaged potential from the Polyakov loop and the anti-Polyakov 
loop at the same z plane. In Fig [5j we plot the resultant potential at T ~ 
167 MeV as for several values of z. We also show the best-fit function of 
— Aexp(— MR)/ R + V . At z = 0, the potential is strongly screened and the 
screening mass is M ~ 1 GeV. As z increases, the screening mass decreases as 
M — > 0. In this phase, we can change the potential from the Debye screened 
form to the Coulomb plus confining form by going away from the membrane, 
without changing temperature. 

In summary, we discussed the four- dimensional Yang-Mills system with a 
three-dimensional fermion membrane, and found that the fermion membrane 
induces a deconfinement layer, which is about 1 fm thick, in confinement 
phase. If we introduce fermions with the larger flavor number, we would ob- 
tain a conformal layer at zero temperature. In general, we can patch not only 



6 



2 




z = 3a 
z = 2a 
z = 1a 
z= 



-A- 

-e- 

■B- 











2 



3 



4 



5 



R [a] 



Fig. 5. The color averaged potential V(R) with N s = 16 and Nt = 6. The curves 
are the best-fit functions of —Aexp(—MR)/R + Vq. 

a membrane but also various different vacuums in the same manner. This kind 
of theoretical study is one possible approach for spatially-nonuniform vacuum 
in lattice gauge theory. 
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